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Abstract. High resolution cosmological N-body simulations show that 
the density profiles of dark matter halos in hierarchical cosmogonies are 
universal, with low mass halos typically denser than more massive ones. 
This mass-density correlation is interpreted as reflecting the earlier forma- 
tion of less massive objects. We investigate this hypothesis in the light of 
formation times defined as the epoch at which halos experience their last 
major merger. We find that the characteristic density and the scale radius 
of halos are essentially proportional, respectively, to the critical density 
of the universe and the virial radius at the time of their formation. These 
two relations are consistent with the following simple evolutionary pic- 
ture. Violent relaxation caused by major mergers rearrange the structure 
of halos leading to a universal dimensionless density profile. Between ma- 
jor mergers, halos gradually grow through the accretion of surrounding 
layers by keeping the central part steady and only expanding their virial 
radius as the critical density of the universe diminishes. 



1. Introduction 



Despite the big efforts made in the last two decades in trying to understand the 
origin of the density profile of galaxies and clusters, the situation is still very 
confused. Is the typical density profile of halos mainly the result of accretion 
(spherical secondary infall) or of repetitive major mergers (strong violent relax- 
ation)? Does it arise, on the contrary, from secular evolution due to the effects 
of dynamical friction plus the tidal stripping of captured satellites? Or is it some 
combination of the preceding processes? 
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Recently, two empirical facts have emerged from high resolution cosmologi- 
cal N-body simulations which seem to have important implications on this long- 
standing debate. First, the spherically averaged density profile of dark-matter 
halos is universal (Navarro, Frenk, & White 1997, hereafter NFW; Cole & Lacey 
1997; Tormen, Bouchet, & White 1997), its exact form being however contro- 
versial (see Moore et al. 1998). The various laws proposed can be encapsulated 
within the general expression 

p(£) & m 

Pcrit e*(l + ^)(7-«)//3 ' U 

with a = 1, /3 = 1, and 7 = 3 for the NFW profile, and a = 1.4, (3 = 1.4, 
and 7 = 2.8 for the Moore et al. (1998) profile. The classical King law and 
the Hernquist (1990) profile also admit this general form for a = 0, j3 = 2, and 
7 = 3, and a = 1, /3 = 1, and 7 = 4, respectively. In equation (||), £ = r/r s is the 
radial distance to the halo center in units of the scale radius r s , and S c = p c /Pcrit 
is the characteristic halo density in units of the critical density of the universe. 
The dimensional parameters p c and r s are linked by the condition that the mean 
density within the virial radius R of a halo of a given mass is a constant factor 
a times Pent! here we adopt a = 200. Therefore, the density profiles of halos 
at a given epoch depend on their mass M through a unique parameter, 5 C or 
x s = r s /R. Second, the smaller the mass of halos, the denser they are (NFW, 
Cole & Lacey 1997) or, equivalently, 5 C is a decreasing function of M. 

This mass-density correlation is interpreted as reflecting the fact that, in 
hierarchical clustering, less massive halos form typically earlier when the mean 
density of the universe is higher. According to this interpretation the structural 
properties of halos would be fixed at their time of formation. As major mergers 
yield a substantial rearrangement of the system, while accretion only causes its 
smooth evolution, this would point to the fact that violent relaxation taking 
place after major mergers is the main responsible for the halo density profiles. 
NFW have shown, indeed, that the characteristic density of halos with a given 
mass at a given epoch is proportional to the mean cosmic density when they form. 
This connection between the characteristic density of halos and their formation 
time strongly favors the previous interpretation of the mass-density correlation. 
However, the physics behind such a proportionality is unclear. Moreover, the 
halo formation time adopted by NFW does not correspond to the last major 
merger experienced by these systems. 

Here we report the results of a similar analysis of the mass-density corre- 
lation recently performed by Salvador-Sole, Solanes, & Manrique (1998, SSM) 
using a better suited formation time estimate. Not only does our analysis con- 
firm the origin proposed by NFW for the mass-density correlation, but also pro- 
vides a physical base for the proportionality between halo characteristic density 
and cosmic density at formation. Furthermore, our results point to a complete 
picture for the evolution of halo structure as a consequence of the interplay of 
accretion and major mergers. 
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2. A Better Suited Halo Formation Time Estimate 



To follow the formation and evolution of halos, SSM have used a slightly modified 
version of the extended Press-Schechter (1974, PS) clustering model, known 
to agree with N-body simulations. The only change introduced consists on 
including a schematic distinction between minor and major mergers. Such a 
distinction does not obviously alter the success of the basic model while it allows 
one to define the formation of a halo as the last major merger it experiences. 

We say that a halo of mass M experiences a major merger and is destroyed 
when the relative mass captured by a halo AM/M, with AM = M' — M the 
increment of mass, exceeds a certain threshold A m . Otherwise, the event is 
regarded as an accretion and the capturing halo keeps its identity. Consequently, 
the mass-accretion and destruction rates of halos of mass M at t are respectively 
given by 

/■M(A m +l) 

rf nass (M, t)= AM r£fe(M -> M', t) dM' (2) 

Jm 



and 



poo 

r d (M,t)= r£ c (M^M',t)dM', (3) 

JM(A m +l) 

where r™ c (M — > M',t) is the specific merger rate in the usual extended PS 
model (Lacey & Cole 1993). 

The distinction between minor and major mergers does not modify the mass 
function of halos at t, N(M,t), which therefore takes the standard PS form. 

As a result of major mergers, some halos are destroyed and other halos 
form. The formation rate of halos with mass M = M(t) can then be obtained 
from the conservation equation 

r i[M(t),t] = +AM(t),t] +d M r* mass (M,t)\ M=M{t) (4) 

for the number density of halos per unit mass along mean accretion tracks, M(t), 
solution of the differential equation dM/dt = r^ ass [M(t), t\. 

Finally, from this formation rate, one can compute the distribution of for- 
mation times for halos of mass Mq at to > t (see SSM for details) 

$ f (t) =r i [M(t),t] expj-^V[M(0,t']dt'j, ( 5 ) 

with M{t) the mean accretion track satisfying M(to) = Mq. The median 
value for this distribution function is then taken as the typical formation time 
tf(Mo,to) of such halos. 

At this stage, the value of A m establishing in a simple way the frontier 
between accretion and major mergers must be regarded as a phenomenological 
parameter of the model. Its value will be fixed in the next section through fits 
to the empirical mass-density (or mass-radius) correlation, the only data which 
seems to be connected with the distinction between merger and accretion. 
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Figure 1. Predicted 5 c (Mq) correlations compared with the empirical 
data from NFW's N-body simulations (filled circles). Dotted and short- 
dashed curves show the predictions for two extreme values of A m , while 
the solid curves correspond to the value of this parameter that gives 
the best overall fit. Cosmogonies with Qq < 1 contain a fourth long- 
dashed curve which shows, for A m = 0.6, the predictions arising from 
the assumption that p c is proportional to p C rit[z{ (-^o)], instead of to 
the mean density of the universe at halo formation. 
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3. The Mass-Density and Mass-Radius Correlations 



As shown in Figure |], the empirical 5 c (Mq) correlation for present halos is well 
fitted, in all the cosmogonies analyzed, by the simple proportionality between 
the halo characteristic density p c and the mean density p(zf) of the universe 
when halos form proposed by NFW 

5 C = Cn [l + Zi (M )}\ (6) 

provided that the mass threshold for merger takes the very natural value of 0.6. 

Similarly good fits are also obtained for this same value of A m in low £1 
cosmogonies if p c is taken instead proportional to the critical density for closure 
Pcir(^f) at the time of halo formation. This new model, which assumes that the 
dimensionless characteristic density of halos at their time of formation is equal 
to the cosmogony-dependent constant 5 c f, while the value of its dimensional 
counterpart p c remains fixed since that epoch, is described by the relation 

In fact, the similarity between the empirical distribution function of halo 
characteristic densities and the theoretical one implied by this latter relation, 
suggests that eq. (f^) holds also for individual halos. More importantly, in con- 
trast with the relation (^) proposed by NFW, eq. (^) has a clear physical inter- 
pretation: violent relaxation occurring after major mergers produces virialized 
halos with identical dimensionless density profiles, which remain essentially un- 
altered during the accretion phase, except for the continuous stretching of their 
outer boundaries caused by the secular decrement of the critical density of the 
universe. 

Given the tight relation between the characteristic density and the scale 
radius of halos of a given mass at a given epoch, if the evolutionary scheme 
we are proposing is correct, the value at formation of the dimensionless scale 
radius x s of halos should also be universal, while the associated dimensional 
parameter r s should keep its value invariant between major mergers. Under 
these circumstances, the scale radius r s for current halos of mass Mo should be 
proportional to the virial radius of the halos when they formed, which implies 
that the empirical mass-radius correlation associated with the mass-density one 
should be well fitted by the expression 

_ R{M[z { (M ),M }} 
Xs ~ Xsl R(M ) ■ W 

We want to stress that this is a new relation, independent from the two previous 
ones (eqs. || and J?]]). Indeed, to verify it we need to estimate the mass of 
halos at the time of their formation (i.e., when they experienced their last major 
merger), so a clustering model such as the one developed in the preceding section 
is necessary. In contrast, the verification of the two previous relations does not 
require any clustering model. 

As can be seen from Figure ^ relation f^) gives very good fits, indeed, to 
the empirical correlations for all cosmogonies, provided, once again, that A m is 
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Figure 2. Predicted x s (Mq) correlations compared with the empirical 
data from NFW's N-body simulations (solid circles). 
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equal to 0.6. Furthermore, as shown in Table |], the best fitting values of x B f are 
in very good agreement with the values of this parameter implied independently 
by the best fits obtained from relation (0) above. These results therefore give 
strong support our physical interpretation of the fitting formulae. 



Table 1. Dimensionless parameters in various Cosmogonies 



P(k) 


n 


Ao 


0-8 


c 


<5 c f 


x si a 




SCDM 


1.0 


0.0 


0.63 


1.21xl0 4 


1.21xl0 4 


0.173 


0.229 


ACDM 


0.25 


0.75 


1.3 


4.21xl0 3 


3.77xl0 3 


0.291 


0.285 


n = —1.5 


1.0 


0.0 


1.0 


8.30 xlO 3 


8.30xl0 3 


0.204 


0.223 


n = -1.0 


1.0 


0.0 


1.0 


1.28xl0 4 


1.28xl0 4 


0.169 


0.181 




0.1 


0.0 


1.0 


2.65xl0 4 


l.OOxlO 4 


0.188 


0.184 


n = -0.5 


1.0 


0.0 


1.0 


2.16xl0 4 


2.16xl0 4 


0.135 


0.148 


n = 0.0 


1.0 


0.0 


1.0 


6.19xl0 4 


6.19xl0 4 


0.088 


0.096 




0.1 


0.0 


1.0 


5.77xl0 5 


1.33xl0 5 


0.064 


0.065 



"implied by S c { 

6 from direct fits to the empirical mass-radius relation 



4. Summary and Discussion 

For A m ~ 0.6, the empirical mass-density obtained in a number of different 
cosmogonies is very well fitted by the simple relations (|6|) or (^). This confirms, 
with a physically motivated halo formation time estimate, the previous claim 
by NFW that the characteristic density shown by dark halos in equilibrium is 
proportional to the mean or critical density of the universe at the time they 
form. A relation like eq. (0) is naturally expected in the case that halos show a 
universal dimensionless density profile at the time of formation, and then keep 
the central part steady while the virial radius expands as the density of the 
universe (and its critical value) diminishes. This simple interpretation of the 
mass-density correlation is supported by the fact that relation (|8|) also gives, for 
A m ~ 0.6, very good fits to the empirical mass-radius correlation, with values 
of x s f fully consistent with those of <5 c f found in the previous fit. 

Such a behavior for the scaling parameters of halo density profiles is very 
well understood as the combined result of the different dynamical effects of ma- 
jor mergers and accretion. The strong violent relaxation due to the former would 
produce an important rearrangement of the halo structure leading to essentially 
a universal dimensionless density profile only dependent on the particular cos- 
mogony considered. On the contrary, the smooth growth via spherical secondary 
infall due to accretion would keep the halo central body unaltered and gradually 
build an extended envelope following the slow decrease of the mean density of 
the universe. 

In another talk in this meeting (see Gonzalez-Casado, Raig, & Salvador-Sole 
1998) we report the results of the detailed check of such an evolutionary scheme 
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for halos. As it is shown there, the values of parameters 5 C { or x s f are completely 
fixed in such a scenario, and it turns out that they are consistent with the values 
found (see Table [j]) in the fits to the empirical mass-density and/or mass-radius 
correlations. 
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